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Abstract 

Predictions for semi-leptonic decay rates of the r lepton into two meson final 
states K~TT iy T , K°tt~v t , K°K~v t , and three meson final states K~ir~ K + u T , 
K°tt-K^u t , K s k-K s v t , K s n-K L v T , K L ir~K L u T , K--k°K°u t , tt it K-v t , 
K~Tt~Tt + v T , 7r~ K°ir°v T are derived. The hadronic matrix elements are ex- 
pressed in terms of form factors, which can be predicted by chiral Lagrangians 
supplemented by informations about all possible low-lying resonances in the 
different channels. Isospin symmetry relations among the different final states 
are carefully taken into account. The calculated branching ratios are com- 
pared with measured decay rates where data are available. 

PACS numbers: 13.35,+s, 11.30.Rd, 12.40.Vv 
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*We wish to dedicate this work to Roger Decker whose death at an early age is a great loss to his 
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I. INTRODUCTION 



With the experimental progress in r decays an ideal tool for studying strong interaction 
physics has been developed. In particular, final states with kaons provide a powerful probe 
of the strange sector of the weak charged current. 

The decay rate for the simplest decay mode with one kaon r — > Kv T is well predicted 
by the the kaon decay constant fx, which can be measured in T(K — > /iz^J. Radiative 
corrections to this decay mode are also available [|IJ. The decay rate for r — > K-nv T on the 
other hand is the Cabibbo suppressed analogy to r — > p{— > 7r7r)z/ T , z.e. the i^7r system is 
expected to be dominated by the if*(892) resonance 0. However, the results for this decay 
rate are fairly sensitive towards the details of the parameterization of the K* resonance, 
such as contributions from higher radial excitations, and we will discuss this dependence in 
some detail in Sec. |H[ Since a reliable parameterization of the K* propagator is also needed 
for the three meson decay modes discussed below, the informations from the r — > K-rv t 
decay rates are used to put constraints on the K* resonance parameterization. 

The hadronic matrix elements for three meson final states with pions and kaons have a 
much richer structure. A general parameterization of the form factors in these decay modes 
was proposed in ||. The physical idea behind this model can be summarized to: 

• In the chiral limit the form factors are normalized to the SU(3)l x SU(3)r chiral 
model. 

• Meson-vertices are independent of momentum. 

• The full momentum dependence is given by normalized Breit-Wigner propagators of 
the resonances occurring in the different channels. Resonances occur either in Q 2 [the 
total invariant hadronic mass squared] which are the three body resonances or in the 
Dalitz plot variables Sj which are two body resonances. 

Parameterizations of the amplitude for the 37r final states within this model can be found 
in [f§-f|. In this case the vector form factor is absent due to the G parity of the pions. The 
pion decay modes offer a unique tool for the study of p, p resonance parameters in different 
hadronic environments, competing well with low energy e + e~ colliders with energies in the 
region below 1.7 GeV. The decay modes involving pions and kaons allow for axial and 
vector current contributions at the same time 0,0]. The vector form factor is related to the 
Wess-Zumino anomaly |||| whereas the axial-vector form factors are predicted by chiral 
Lagrangians as mentioned before. A parameterization of the r~ — > K~tt~ K + v t decay mode 
within this framework has been developed in |K|. However, our result for this decay mode 
differs from the their result (see Sec. |V11| ). 

In the present paper, we reexamine the model used in || for the r decay modes involving 
pions and kaons. We take into account possible isospin and SU(3) symmetry relations, which 
allow for additional resonance contributions to the form factors. We reanalyze the issue of 
the strange axial resonances {K\ states) in view of new experimental results, and we include 
uo — $ mixing. The new parameterization leads to sizable differences in the predictions of the 
decay rates compared to ||. In addition, we derive a parameterization for the final states 
with two neutral kaons r _ — > Ks^~ Ksf T i r ~ ~^ Kltt~ Klu t1 and r~ — > Kg-K' K^v T . 
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The paper is organized as follows: In Sec. II we discuss the two meson decay modes 
A"~7r z/ T , K^tx~v t and K°K~v T and fix the parameters of the K* resonance. In Sec. Ill we 
review the general structure of the matrix elements of the weak hadronic current in the three 
meson case. In Sec. IV (Sec. V) we discuss in detail the individual matrix elements of the 
axial- vector (vector) current, carefully exploiting isospin and flavour symmetry. Final states 
with two neutral kaons are considered in Sec. VI. We give our numerical results in Sec. VII, 
and finally in Sec. VIII we give a brief summary. 



II. THE TWO MESON DECAY MODES 

As mentioned before, a reliable parameterization of the K* propagator is needed for the 
three meson decay modes with kaons. Since the decay modes r — > K~tt v t and r — > K°7r~v T 
are dominated by the K* resonance, we use the experimental informations on the decay rates 
into these final states to fix the parameters of the K* resonance. 

The matrix element for the semi-leptonic decay into two mesons h a and h b 

r(l, s) -> i/(7', s') + h a (qi, mi) + h b (q 2 , m 2 ) , (1) 
can be expressed in terms of a leptonic (M M ) and a hadronic vector current (J M ) as 

M = J= Ct) . (2) 

In Eq. (H), G denotes the Fermi-coupling constant and 9 C is the Cabibbo angle. The leptonic 
and hadronic currents are given by 

M m = «(Z', S ')7m(1-75MM), (3) 

and 

J*{<h,q2) = (ha{qi)h b (q 3 )\V'>{0)\0). (4) 

The hadronic matrix elements for the Cabibbo suppressed decay modes K~ 

are dominated by the K* resonance Tj^}(Q 2 ) } whereas the one for the Cabibbo allowed mode 

K°K~ is dominated by the high energy tail of the p. One has 

< ^-( gi )7r°(g 2 )|n0 > = 77f < ^(?iK(<? 2 )ino > , 

= ^=T^){ qi -q 2 ) v T^ (5) 

< K { qi )K-{q 2 )\V^\0 > = T«( gi - q 2 ) u T^ , 

where Q = q\ + q 2 and T^ u denotes the transverse projector, defined by 

G^G V 
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The form and the normalization of the hadronic matrix elements in Eq. (||) are fixed by chiral 
symmetry constraints, which determines the matrix elements in the limit of soft meson 
momenta Q 2 — > 0. The strong interaction effects beyond the low energy limit are taken 
into account by the vector resonance factors T^\Q 2 ) and T^}(Q 2 ) with the requirement 



Tx\Q 2 = 0) = 1 (X = p, K*). Note that in the case of the K* we have neglected a 
scalar contribution proportional to This scalar part is proportional to the off-shellness 
(m 2 K ^—Q 2 ) of the K* and therefore strongly suppressed. We have checked its size numerically 
and found it to be negligible. 

Note that our results for the hadronic matrix elements for K~n° and K°7r~ in Eq. ([5]) 
differ from the results in the Tauola Monte Carlo by an overall factor. In fact we believe 
that Tauola is off from the correct normalization by a factor of 2/\/3, i.e. the corresponding 
matrix elements in Tauola should be multiplied my v^3/2. Note furthermore that Tauola 
does not include the i\"*(1410) in the Kir resonance, i.e. it corresponds to (3k* = in 
Eq. ([K]). These two differences with our parameterization almost cancel each other, such 
that Tauola (with fix* = 0) gives a number for the r — > K*u T branching ratio which is very 
close to our result with /3 K * = —0.135 [see below]. 

If not stated otherwise, we use the following form for the two particle Breit-Wigner 
propagators with an energy dependent width throughout this paper: 



BW x [s] 



\M 2 x -s-iJST x (s)} ' 



where X stands for the various resonances of the two meson channels. For a 1 
the energy dependent width is H] 



(7) 

2 decay, 



rvM 



r x —(—) 



2n+l 



Px/ 



p = -!W( a - (Mi + M 2 ) 2 )(s - (M 1 - M 2 f) 

L\/ S 



Px 



2M 



lr^J(M 2 x - (Mx + M 2 ) 2 )(M1 - (Mx - M 2 ) 2 ) , 
vi x 



where n is the power of \p\ in the matrix element, i.e. n = 1 for the decay modes in this 
paper. 

We use the following parameterizations for the p resonance: 



BW J s)+(3 p EW As] 



(8) 



where 



0p 



m. 



-0.145 , 
0.773 GeV,r p = 
1.370 GeV,I> 



0.145 GeV, 
= 0.510 GeV. 



(9) 



These are the values which have been determined from e + e — > 7r + 7r in [|J] and have been 
used in |Q for the non-strange case. For the vector resonances with strangeness, only the 
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7T*(892) was considered in ||. Our parameterization for T^}(Q 2 ) allows for a contribution 
of the first excitation K* (1410) in analogy to Eq. (0): 



1+0 



K* 



BW^( S ) + ^BW^(s) , (10) 



where 



m K * = 0.892 GeV , T K * = 0.050GeV , 

m K *> = 1.412 GeV , T K *> = 0.227 GeV . (11) 

In the limit of 577(3) flavour symmetry, one would expect a contribution to the K* resonance 
from the first excitation K* (1410) with the same relative strength (3k* = —0.145 as mea- 
sured in the non-strange case. The numerical results for the decay rates K~ 
are very sensitive to the parameter (3k* in Eq. fllDp . Since a reliable parameterization of the 
K* propagator will be needed for the decay modes considered in this paper, we will use the 
decay mode r — > K*v T to fix the parameter (3k*- 
Adding up the two charge modes in Eq. <ffi) 



we obtain 



B(K*v 


r) = B(K~ 


ttV) 


+ B{K°tt- 




B(K*u T ) = 


1.00% 


for 


Pk* = 


0., 


B(K*u T ) = 


1.28% 


for 


Pk* = 


-0.11 , 


B(K*v T ) = 


1.36% 


for 


13k* = 


-0.135 


B(K*v T ) = 


1.44% 


for 


Pk* = 


-0.157 



(12) 



(13) 

The main effect is due to the normalization factor 1/(1 + (3k*) 2 multiplying the TT*(892) 
contribution, whereas the TT*(1410) contribution is strongly phase space suppressed. 

The latest experimental result on the branching fraction B{K*v T ) is 1.36 ± 0.08 % [12 



Thus our results favour a negative value of (3k*, and we obtain (3k* = —0.135 ± 0.025 as a 
result of this analysis. This value is remarkably close to the strength of the p' contribution to 
the p Breit-Wigner in Eqs. (HH), supporting the use of approximate SU (3) flavour symmetry. 

Note that we use here the values of [|T3] for the mass and widths parameters of the K* 
in Eqs. (|H],|ll]), whereas the p' mass and width parameters in Eq. (||) have been determined 
from a fit to the Q 2 distribution in r — >• 2ixv r [|J. Indeed, a more reliable determination of 
the parameters of the off-shell K* propagator (including (3k*) could be obtained from a fit 
to the Q 2 distribution in r — > K*u T . Unfortunately, there is no experimental information of 
the Q 2 dependence in this decay mode available, and we use therefore the simpler approach 
of fixing (3k* from the branching ratio and using the K*' parameters from |13[ . 

For the decay into two kaons we obtain from the matrix element in Eq. (El) 

B(K°K~v T ) = 0.11% , (14) 



in good agreement with the recent world average B(K K v T ) = 0.13 ± 0.04% |L2 . 
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III. GENERAL STRUCTURE OF THE WEAK MATRIX ELEMENTS IN THE 

THREE MESON MODES 



Let us briefly recapitulate the general structure of the semi-leptonic decay 

t(1, s) -> s') + h a {q 1 ,m 1 ) + h b (q 2 , m 2 ) + h c (q 3 , m 3 ) , (15) 

as introduced in [Q. In Eq. (|H>D, h(qi,mi) are pseudoscalar mesons and in our case of 
interest, at least one of the mesons is a kaon. The matrix element is 

M = ^{Zk)M»J\ (16) 

where the cosine and the sine of the Cabibbo angle {6c) in Eq. ([IB] ) have to be used for 
Cabibbo allowed AS = [i.e. final states with two kaons] and Cabibbo suppressed (AS*! = 1 
[i.e. final states with one kaon] decays, respectively. The leptonic (M M ) current is already 
given in Eq. (|3|) and the hadronic (J M ) currents can be written as 

(72, ft) = (h a (qi)h b (q 2 )h c (q 3 )\V^0) - A"(0)|0) . (17) 

and are the vector and axial-vector quark currents, respectively. The most general 
ansatz for the matrix element of the quark current J M in Eq. ([17]) is characterized by four 
form factors 

^(qu Q2, ft) = VI 1 Ft + Vf F 2 + i Vg F 3 + Vf F 4 , (18) 

with 

Vf =( qi -q 3 ) u T^, 
V$ = (q 2 -q 3 ) u T^ , 

2 W2 q 3 Ju 

V* =e^q la q 2p q 37 , 



A' 



denotes the transverse projector, defined in Eq. (^). Since the strong interaction con- 
serves parity, the axial- vector current induces the form factors Fi, F 2 and F4 while the vector 
current induces F 3 . In the limit of vanishing quark masses, the weak axial- vector current 
is conserved and this implies that the scalar form factor F4 vanishes. The massive pseu- 
doscalars give a contribution to F4, however, the effect is very small |14[ and we will neglect 



this contribution in the subsequent discussion, i. e. we set F 4 equal to zero. All form factors 
Fi are in general functions of Q 2 , Si = (q 2 + g 3 ) 2 , s 2 = (qi + g 3 ) 2 and s 3 = (qi + q 2 ) 2 . 

A specific model for the form factors Fi for various three meson final states was derived 
in [Q. This model takes into account the chiral symmetry constraints of QCD as well as 
the resonance phenomena present in r decays and has been developed in . Let us briefly 
summarize the results (for more details see H]): The chiral symmetry constraints lead to the 
following expression for the hadronic axial- vector current in the limit of soft meson momenta: 

(h a h b h c \A"\0) = 2V lf bC) {GitUqi - q,)u + G^ ft (q 2 - q 3 ) v } , (20) 
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where the coefficients A^ abc ' are given in Tab. | for the various decay modes and the coeffi- 

-,(abc) 
r i,soft 



cients G^l t are all equal to one 



Gg& = l, («• = !, 2) (21) 

except for G^ so J t K ■* and G^ so f t n ■* which vanish: G 2 K so ^ t K ■* = G^ 2s ^ t n ^ = . The vector 
current arises from the Wess-Zumino Lagrangian for the axial anomaly ||||. One obtains 
in the low energy limit: 

(h a h b h c \V»\0) = n ^^ We^^^^S • (22) 



2V2tt 2 /!" 

„™ t^u ITU n( ah 

oft 



The coefficients A^ 6 ^ for the vector current are given in Tab. [TJ and G^^ t = 1 , except for 

G ioft K ° ] and G t°oft K } which vanish: G { £ m fl K ° ) = G { ^ s f ft K ] = . 

As mentioned before, the strong interaction effects beyond the low energy limit are taken 
into account by inserting resonance form factors Gi <2j3 (Q 2 , s») [i = 1, 2, 3] into the amplitudes 
of Eqs. (|20| , p2D with the requirement G^l ~^ G ^2%soft = 1 i n the limit Q 2 , Sj — > 0, except 
for G 2 K 3 n K \ G^ K " ^ and ^ K \ which vanish in the chiral limit. 

In fact, the functions (jx23 are products of normalized Breit-Wigner resonances in Q 2 
and Sj. Comparing Eqs. fl2"0l , |22"D with the general expansion in Eqs. ( |l8l , |l9l) leads to 

Fi abc) (Q 2 ,s 2 ,s 3 ) = 2 ^f bC) G[ abc \Q 2 ,s 2 ,s 3 ) , (23) 

3 J TV 

F 2 iabc \Q 2 , Sl ,s 3 ) = 2 -^-^ G 2 abc \Q 2 ,s l ,s 3 ) , (24) 
Ft bc \Q\ 8l , s 2 , s 3 ) = ^^j- 3 G { 3 abc \Q 2 , *i, « 2 , s 3 ) • (25) 

The Breit-Wigner functions G\ j2 (G 3 ) are listed in Tab. [I] (|T|) for the various decay modes 
and the precise form will be discussed in the subsequent sections. Note that by convenient 
ordering of the mesons, the two body resonances in F\ (F 2 ) occur only in the variables s 2 , s 3 
(si,s 3 ). 



IV. ISOSPIN AND FLAVOUR SYMMETRY RELATIONS IN THE 
AXIAL- VECTOR CURRENT 

Amongst the hadronic matrix elements for the hadronic final states K~ir~K + , K°n~K° 
and K~n°K° (which are versions of K71K with different third components of the isospin) 
and amongst those for 7r 7i K^ , K~ty~'k + and n~K 7r° (different versions of Kttit), there are 
isospin symmetry relations which have not been fully exploited in |y. Consider for example 
the final state K~tt~K + , which is dominated by a (K*) resonance in ir~K + . In Tab. [I] of 
@, this is taken into account by 

g (k-,-k+) = BWai (Q 2 )T^( Si ) . (26) 
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However, by isospin symmetry, the amplitudes (K*)° — ► tc°K° and (K*)° — > 7r fT + have a 

ratio of Jl/S : (— )y / 2/3. Taking into account the different normalization coefficients A^ abc \ 
this immediately leads to 



bw Ai (q 2 )[It k *( Si ) + -- 



(27) 



which differs from 71 K "* = in Tab. [I] of p3]. The dots indicate that there must be an 
additional contribution which ensures that ^ K — > for Q 2 , Sj — > as required by the 
chiral limit (see below). 

In fact, by using all the relevant Clebsch-Gordon coefficients, we can predict the matrix 
element for K~tc°K° from those for K~ix~K + and K°n~K°, and the one for 7i~K°ii can 
be predicted from those of K and K n z . In using the Clebsch- Gordons, some care 
is needed because in ||, the phase conventions are not identical to the Condon- Shortley-de 
Swart ones. And so, in a first step Tab. [I] of has to be transformed into Condon-Shortley- 
de Swart by multiplying A (abc) with (— 1) for each tt~ , K~ , K°, K° in the final state. Then 
the Clebsch-Gordans may be used, and finally the results are transformed back into the 
conventions of @. This issue does not affect the above example, because < K~ir~K + \ and 
are even under this transformation. 

Taking the matrix elements for K~n~K + and K°n~K°, K and K n n + from ||, 
we obtain (see also Tab. ITf) 



x 



3/ w 2x/2 



(28) 



< ir~ K°7r°\A>*\0 >- 



2^3 sin £ 

3fn 2^2 



BW Kl (Q 2 ) T» v 



X 



h - q3)u + 



Tic*(si) - Tjr*(s 3 ) (g 2 - q-s 



(29) 



Note that the form factor G2 in the two decay modes K~n K and n~K ir receive con- 
tributions from Tk* resonances, which however vanish in the low energy chiral limit, where 
our Breit-Wigner propagators are normalized to one. In @, this contribution was not taken 
into account. 

The parameterizations of the two particle resonances have been determined in Sec. ||, 
where we found (3k* = —0.135 ± 0.025 for the K* resonance. Our numerical results for the 
three meson decay modes are very sensitive to the parameter (3k* and we will discuss this 
dependence in more detail in Sec. |VII| . 

Let us now discuss the three particle resonances. As in ||, we use the A\ resonance in 
the non-strange case with energy dependent width 



BW Al (s) 



ni 



im Al T Al g(s)/g(m Al ) 



(30) 



8 



m Al = 1.251 GeV , T Al = 0.475 GeV . (31) 

where the function g(s) has been calculated in j3|. Note that we use a smaller A\ width 
than used in |3|J5|,§] (T Al = 0.599 GeV). The formalism of structure functions, which was 
developed in Hfl, allow for a much more detailed test of the hadronic matrix elements 
and the resonance parameters than it is possible by a rate measurement alone. However, 
the ratios of structure functions, which were predicted in [^]|| for the r _ — > it~it~-n + v T 
decay mode and found to be in good agreement with the experimental data JIB| , are not 



sensitive to the A\ mass and width parameters since they cancel in the ratios. On the other 
hand, a measurement of the Q 2 dependence of the structure functions itself would be very 
sensitive to the A\ parameters jjT7[] (after having fixed the p parameters in the ratios). An A\ 
width of 0.475 GeV predicts a branching fraction B(r~ — > tx~tx~tx + v t ) of 8.6 % (compared 
to 6.3 % for F Al = 599 GeV), which seems to be favored by the experimental results for 
B(3hu T ) = B(n-ir-ir + is T )+B(K-n-n + is T )+B(K-ir-K + v T ) + B(K-K-K + is T ) = 9.51 ±0.18 

0- 

There are two relevant particles which have the right quantum numbers in the case of 
the three particle resonances with strangeness, i.e. the i^i(1270) and the ^(1400). Results 
from the TPC/Two-Gamma Collaboration had indicated that only the i£i(1400) occurs in r 
decays into Ktttt [JTSfl . Thus in || only the i^i(1400) was included. Recent experimental data, 
however, do not support the results from ||18j| , but rather indicate that both K\ resonances 
are produced, with substantially more i^i(1270) than ^(1400) ||19|| . In fact, this is much 
more natural within the framework of the resonance enhanced chiral calculation. There are 
two types of decay chains which lead to the Ktttt states, i.e. 

a) K x K*tt, K* -> Kir , 

b) K\ — > pK, p — >• mi . 

The i^i(1400) decays exclusively into K*7T, while the Xi(1270) decays mainly into Kp, with 
a smaller branching ratio into K*tt. Thus one should assume that the decay chain (a) is 
dominated by the ^(1400) with a little admixture of i^i(1270), whereas in (b) only the 
7^(1270) contributes: 



BW Xl (i400)(s) +£BW Kl (1270)(s) 



1 + e 

TP i (s)=BW Kl{1270) (s). (32) 

Using the ratio T(Ki(1270) ir*7r)/r(ir 1 (1400) -> K*7r), we obtain after phase space 
correction: 



If I = 0.33. (33) 

Thus £ is determined up to a sign ambiguity. We will consider both possibilities below and 
find that £ = —0.33 leads to decay rates which strongly disagree with experimental data 
(see Sec. |VII|) . Therefore we believe that £ = +0.33 is the physical choice. In case of the K\ 
resonances we use normalized Breit-Wigner propagators with constant widths: 
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BW Kl [s] = K l , Kl 34 

[s -m 2 Ki +im Ki r Kl ] 



and 13 



m Kl {UOO) = 1.402 GeV, 1^(1400) = 0.174GeV, , , 

m Xl (1270) = 1.270 GeV , T Xl (1270) = 0.090GeV . ^ > 

All matrix elements of the axial weak current are summarized in Tab. |L The parameteriza- 
tion for the decay modes with two neutral kaons will be discussed in Sec. [VI]. 



V. ISOSPIN AND FLAVOUR SYMMETRY RELATIONS IN THE VECTOR 

CURRENT 

We consider now the vector current with its contributions from the anomaly. Let us first 
correct an error in Tab. |I]of ||. The two kaon resonance in the decay channels K~-k~K + 
and K°Ti~W can not be a p resonance. Due to G parity conservation, the vector resonance 
V in the decay chain W~ -»• p~, p" -> Vn', V -> (K~K + , K°K°) must have G = -1. Thus 
only the to and the $ qualify. In the limit of exact SU (3) flavour symmetry, the coupling 
p$7r vanishes, such that only the uj contributes in this limit. Taking u;$ mixing into account, 
we use 

T u (s) = -i-[BW w (s) + eBW^(s)] , (36) 



where 20 



e = 0.05 . (37) 

Obviously, use of the uj instead of the p changes the relevant isospin symmetry relations. 

We use Breit-Wigner propagators with fixed widths, of the same form as the one used 
for the Ki resonances in Eq. (|3"4D, with [T^ 



m w = 0.782 GeV , T w = 0.00843 GeV , , 
= 1.020 GeV , 1^ = 0.00443 GeV . ^ ' 

Similarly to the case of the axial- vector current, there are also different decay chains 
in the vector current case which have either K* or p/w two particle resonances. Whereas 
in the case of the axial- vector current, the relative strength of the strangeness S = and 
S — —1 resonances could be fixed by the chiral limit in 0, an additional free parameter a 
was introduced in the case of the vector current, parameterizing the relative strength of the 
two body resonances. In the case of the K~it~K + and K°tt~K° states, the parameterization 
in was 

G iK~,-K +) = G (K°«-K*) = T ( >2)m n i) (s2)+aT K .(S l ) 

' 1 + a 

where T^ 1 ' should read T w , as we have explained above. However, by writing the relevant 
vertices in a flavour invariant way, 
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gvv^ ua hi(d^V v d a VpIl) 



Wvt, 



T tr(^[^n,n]) (40) 

(where II = 7r a A a describes the pseudoscalar mesons, and similarly V p parameterizes the 
vector mesons) the parameter a can be determined, with the result 

1/V2. 



a 



(41) 



Furthermore, as in the case of the axial- vector current, there are isospin symmetry relations 
between the matrix elements. 

Our final result for the matrix elements of the vector current, taking into account all 
isospin and flavour symmetry relations, can be found in Tab. |J. Note that there are two cases 
in which the anomaly contribution vanishes in the low energy limit, viz. K-7T°K (because 
"anomalies do not develop second class currents", see also 0) and K (because of Bose 
symmetry and the antisymmetry of the anomaly). In both cases, however, the presence of 
resonances leads to non-vanishing contributions in the higher energy regime. Their precise 
forms are predicted from other matrix elements which are non zero in the chiral limit and 
to which they are related by isospin symmetry. These contributions were not taken into 
account in the parameterization used in |§. 

Let us now discuss the three particle vector resonances T p (2) and T$. In [§, a form 
for including p, p' and p" was used, which was obtained from a fit to e + e~ — > rjim 
data plT,|i~0"| . In the three particle vector resonance with strangeness, only the if* (892) was 
included. However, the higher radials will also be included in this paper and we will therefore 
use 



T (2) 
P 



1 



T, 



(2) 



1 + X + p 

1 

1 + X + p 



BW p (s) + \BW p ,(s) + /iBW p »(s) 
+ ABW x *'(s) +pBW K 



BW 



A* 



(s) 



where 



For the p's, 



A 



6.5 
^26 



-0.25 



0.773 GeV 
1.500 GeV. 
1.750 GeV 



1 



-26 



-0.038 



(42) 



(43) 



0.145 GeV, 
: 0.220 GeV, 
= 0.120 GeV. 



(44) 



This is exactly the parameterization of [ID|, which was used in [HI, written in a slightly 
different way. For the K*'s, we use the values ||13||: 



m K * 
m K *> 
m K *" 



0.892GeV,r A -* 
1.412 GeV, T K + 
1.714 GeV, T K * 



0.050GeV , 
0.227 GeV, 
= 0.323 GeV . 



(45) 
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We use energy dependent widths here. 

Similar remarks as have been made above for the T%1 apply here as well. The parameters 
A and p, as well as the p' and p" parameters of have been obtained from a fit to data. 
For the K* in the anomalous channel, we use the same parameters A and p as for the p 
with K* , K*" parameters taken from ]TJ|. In principle, a more reliable determination of 
T^l could be obtained from a fit to suitable data. It should be noted, however, that the 
numerical significance of these details is fairly small, because of the small vector channel 
contribution to the relevant decay modes (see the last 3 entries in column 3 of Tab. PTI| ). 

VI. FINAL STATES WITH TWO NEUTRAL KAONS 

In || the hadronic matrix elements have been expressed in terms of the strong interaction 
eigenstates K° and K°. The actual measurements, however, are performed in terms of the 
weak interaction eigenstates K$ and Kl- Neglecting CP violation, the relation between 
these are given by 

Kq = = — Kj = = — 

K L + K S K L -K S 

Ko = —sr- Ko = —sr- (46) 

Since the weak current J M produces states with |,S| = 0, 1 only, i.e. 

< K°n-K°\J tI \0>=<KVn-KV\J tl \0>=0 (47) 
we can easily show that 

< K L (q 1 )7T-(q 2 )K L (q 3 )\J fi \0 >= - < K s (qi)ir-(q 2 )K s (q 3 )\J»\0 > 

= ^{< K°( qi )7r-(q 2 )KZ(q 3 )\J»\0 > + < J^tT ' (q 2 )K°(q 3 )\J»\0 >} , 

< K s (q 1 )Tr-{< h )K L (q 3 )\J»\0> 

= l -{< K ( qi )7r-(q 2 )KZ(q 3 )\J»\0 >-< (q 2 )K\q 3 )\J^ >} . (48) 

Let us start with the discussion of the axial part of the weak current. From Eq. (GO) and 
Tab. |, we have 

< k\-W\a^ >= ^ Bw^(g 2 ) 

* {Tj; 1] (s 2 )( qi - q 3 ) u + T K *( Sl )(q 2 - q 3 )„} . (49) 
Using Eq. fl4~8"| ) one obtains 
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< K L ( gi )n-(q 2 )K L (q 3 )\A»\0 >= - < K s (qi)K~ {q 2 )K s (qz) K|0 >= 
<K s ( qi )n-(q 2 )K L (q 3 )\A»\0>-- 



cos 6 



2V2 / 
3/ w V 4 



'c 



Bw Al (g 2 



x 



{ [2T«( S2 ) + r«( S3 )] (9i - ft), + [T$(ai) - rW( S3 )] - ft),} ^ . 



(50) 



We turn now to the matrix element of the vector current. According to Tab. |, the relevant 
matrix element is 



< K U 7T 



0_-~ 



-K°\V»\0 >= —^—oosBjf {Q 2 ){V2- l)[V2T w {s 2 ) +T$( Sl )}e^q la q 2f} q 3j . 

(51) 

In this case Eq. ([48]) yields (see Tab. |J) 



< Kltc-KlIV^O >=-< Kgir-KslV^Q > 
i cos 9 C 



2^2/3 2 



:(V2-l)Tf(g 2 )[T^( Sl ) -T K *(s 3 



e^ a(3l q la q 20 q^ 



< K s tt-K l \V"\0 > 
i cos 9 



2^2tt 2 L 3 2 



(V2 - l)Tf (Q 2 ) [2V2T^s 2 ) + T$( Sl ) + 4 1 2( S3 )]e^ 7 giag 2 ^3 7 • (52) 



Whereas the relative amounts of K$Ks and K^K^ states are fixed by general symmetry 
considerations to be equal 



T(K s n-K s ) 
T(K L 7t-K L ) 



(53) 



the ratio R of KqKq and KqKl states 



R 



T(K s n-K & 



T(K S 7i-K L ) ' 

is model dependent and will be discussed in the next section. 



(54) 



VII. NUMERICAL RESULTS 

After having fixed our model for the form factors, we next present numerical results for 
the hadronic decay widths T(abc) normalized to the leptonic width T e and for the branching 
ratios in Tab. |TJ. To calculate the branching ratios, we use the theoretical prediction for 
r e /r^ ^ = 17.8% based on the experimental values for the tau mass m T = 1.7771 GeV and 



lifetime r T = 291.6 fs [22|, rather than using the experimental branching ratio. 
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The decay rate for r decays into three mesons can be calculated from 



12m T Vsm ^ (4vr 



where 



W^4 = (x? + ^)|F 1 | 2 + (x2 + ^)|F 2 | 2 + 2(xix 2 -x2)Re(F 1 F 2 *) (56) 

W B = x 2 |F 3 | 2 (57) 

The variables Xi are defined by X\ = V* — qf — , x 2 = V$ = q% — q%, x 3 = — qf — 

— = V3 2 = \/Q^x 3 q^, where qf (gf) denotes the x (y) component of the momentum 
of meson i in the hadronic rest frame as introduced in They can easily be expressed 

in terms of si, s 2 and S3 0. Eq. (^) shows that there is no interference between the axial 
vector contributions (Fi,F 2 ) and the vector current contribution (F3) in the total decay 
width. 

Our numerical results for the normalized decay widths T(abc)/T e and the branching 
ratios B(abc) for our preferred parameter choices are given in Tab. |T| for the various decay 
channels abc. To get a feeling for the numerical importance of the vector current (i.e. 
the "anomaly"), we list its contribution to the decay width in column 3 of Tab. [TTT]. For 
comparison, we have also listed the available experimental data in column 5 of Tab. |T|. 

Our numerical results for B(Ksir~Kg), B(K~-k~-k + ) and B(7r~ K°7i°) appear to be con- 
siderably higher than the experimental results, whereas the other predictions agree fairly 
well. 

Our results for the K~tt~K + final state do not agree with the results in [Hj, where 
the contribution of the axial-vector channel amounts to less than 10% to the decay rate in 
this channel. In fact, our predictions for the axial-vector contribution is about 60%. This 
result is fairly insensitive towards the details of the K* parameterization (see Tab. p||V| 



and 0, as discussed below). It is however sensitive towards the A\ parameters. Use of 
Ta 1 = 0.599 GeV in Eq. (El|) reduces the axial-vector contribution to about 48%, which is 



still considerably larger than the value in [JTO ] 



Let us now comment on the sign of parameter £ in the parameterization of the strange 
axial resonances K\ in Eq. (P3f), which is only determined up to a sign ambiguity. Our 



predictions for the Cabibbo suppressed decays in Tab. |T| are obtained for £ = +0.33. For 
comparison, Tab. |^ shows the results for £ = —0.33. The choice £ = +0.33 is clearly 
preferred by the experimental data and we therefore believe that this is the physical choice. 

Next we discuss our choice of the two body vector resonance with strangeness, i.e. the 
parameterization of T^l in Eq. (|i~0"p. We have used a X*'(1410) contribution in T^} with 
a strength of (5k* — —0.135 relative to the if* (892), as determined from r — > K*v T (see 
Sec. II). Use of (5k* = —0.11, which is also consistent with the r — > K*u T decay rate leads 



to the results in Tab. 0. The results are fairly close to the numbers presented in Tab. |TT 



Use of (5k* = as in leads to the results in Tab. |VT|, which overall agree better with the 
experimental results than those obtained with (5k* = —0.135. 

At this point a few comments are in order. Firstly, by fixing the coupling constants 
from the chiral limit, we effectively assume exact SU(3) flavour symmetry for the coupling 
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constants. The most important effects of flavour symmetry breaking have been taken into 
account by using the physical masses and decay widths of the mesons in the propagators 
and in the phase space. However, one could expect additional explicit flavour symmetry 
breaking in the couplings. The approximation made may lead to an error of about 10- 
30% on the matrix element level. Thus we consider the agreement between theory and 



experiment in Tab. [HJ as reasonable. Note that this approximation will mainly lead to 
an overall normalization error of the form factors rather than to a modification of their 
momentum dependence in the relevant physical region. 

Secondly, we remind the reader that a change of (3k* has two different effects. It changes 
the normalization of the K*(892) contribution, which is proportional to 1/(1 + 13k*) 2 in the 
rate, and it changes the size of the X*(1410) contribution, which is however strongly phase 
space suppressed. The main effect of choosing a different (3k* is to change the normalization, 
which may well compensate the error being made by deriving the coupling constants using 
flavour symmetry. 

Thirdly, we wish to emphasize that in spite of this normalization uncertainty of the 
matrix elements, our parameterizations are by no means arbitrary. We believe that they 
give a very good description of the resonance substructures which determine the decay rates 
as well as differential distributions and structure functions. 

Our choice for Ya 1 = 0.475 GeV was already discussed before [see the discussion after 
Eq. (0)]. An Ax width of = 0.599 GeV decreases the branching fractions for the K-kK 
decay modes in Tab. [H| by about 15%, which is of course entirely due to a decrease of the 
axial- vector contribution (see Tab. []]). 

As already mentioned in Sec. |VT], the decay rates for Ks7r~Ks and Kltt'Kl are the 
same. This follows immediately from Eqs. (|46|j47|) and is a model independent statement. 
On the other hand, the ratio R 

R ~ T(K s n-K L ) ' (58) 

is model dependent and we obtain from Tab. |TXT| 

R = 0.48 (59) 

This number is fairly insensitive towards a variation of the K* resonance parameters. How- 
ever, use of T Al = .599 GeV in Eq. ([31]) would lead to R = 0.42. 

We would like to point out that a study of angular correlations of the hadronic system 
allows for much more detailed studies of the hadronic charged current (including the details 
of the two and three body resonance parameters) than it is possible by rate measurements 
alone. Of particular interest is the angular distribution of the three mesons in the three 
meson rest frame. The distribution of the normal on the hadronic plane with respect to ul 
(the direction of the laboratory as seen from the hadronic rest frame) allows for a model 
independent separation of the axial- vector and the vector current contribution, /i.e. the 
structure functions Wa(Q 2 , si, S2, S3) and Wb(Q 2 , si, S2, S3) in Eq. ([55]) can be determined 
separately even without reconstructing the r rest frame ||. An experimental analysis of the 
Q 2 and Sj distributions of these structure functions would clearly help to test the parameter- 
izations in Tab. | and Tab. |H| unambiguously. More general distributions like the rotation of 
the mesons around the normal, allow for even more detailed studies of the hadronic matrix 
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element ||. Q 2 distribution for structure functions for K~ir~K + , K~ix~ix + and 7]n^n + final 
states have already been presented in [23] based on the parameterization in ||. We will 
study the Q 2 and the full Dalitz plot distributions for the decay rates and the two structure 
functions Wa and Wb based on the model in this paper in a future publication. 



VIII. SUMMARY AND CONCLUSIONS 

We have discussed tau decays into final states with one or two kaons. The decays into 
ttK are dominated by the K* resonance and therefore allow for a determination of the 
parameters of the K* propagator. The experimental branching ratio can be used to obtain a 
rough estimate of the K*(1410) contribution, but we would like to urge for a detailed study 
of the invariant mass distribution of the hadronic system in order to measure the strength 
of its contribution. The decay into two kaons is predicted in agreement with experimental 
data assuming dominance by the high mass tail of the p. 

The three meson final states KttK and tcttK allow for a much more involved resonance 
substructure. We have extensively reanalyzed these on the basis of the model of |3J. Our 
final results for the branching ratios with our preferred parameter choices have been given 
in Tab. III]. They compare reasonably with experimental data, but are rather sensitive to 
parameters such as (3k*- 

In the case of two neutral kaons, we have expressed the matrix elements in terms of the 
Ks and Ki states and given a prediction for the model dependent ratio of the Ks^~Ks and 
K s ii~K L final states. 
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TABLES 



TABLE I. Parameterization of the form factors F\ and F 2 in Eqs. (23 
elements of the weak axial-vector current for the various channels. 



for the matrix 



channel 

(abc) 



j^(abc) 



G[ abc \Q*,s 2 ,s 3 ) 



G 2 \Q 2 ,Si,S 3 ) 



K--R-K+ 



cos6 r 



Ks-k-Ks 
K s n~K L 
K-n°K° 



—tt — COS 6 r 



2 

— cos C 

4 

— COS C 

4 

3 cos C 
2V2 



BW Al (Q 2 )T^(s 2 ) 
BW Al (Q 2 )T^(s 2 ) 
BW Ai (Q 2 )t£!(s 3 ) 



BW Ai (Q 2 )t£}( Si ) 
BW Al (Q 2 )T$( Sl ) 
-BW Al (Q 2 )[T^( Sl )+T^(s 3 )} 



BW Al (Q 2 )[2TW(s 2 ) +t£}(s 3 )} BW Al (g 2 )[T«( Sl ) -T$(* 8 )] 



BW Al (Q 2 



^(s 2 ) + ^}(s 3 ) 



1 



^BW Al (Q 2 ) r#(.s : ,) 



K 7T 7T 



7T-K 7T C 



sin 6* c 
4 

— sin # c 
2 

3 sin C 
2V2 



T Ki(Q 2 ) T K*( S 2) 



T Ki(Q 2 ) T K*(Sl) 



Tg(Q 2 )T$(s 2 ) 



Z -n]{Q 2 )T^\s 2 ) 

+\t { £{Q 2 )t${s 3 ) 



t kI(Q 2 ) Z#(*i)-T#(*s) 



^(1) 
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TABLE II. Parameterization of the form factor F3 in Eq. ( |25| ) for the matrix elements of the 
weak vector current for the various channels. 



channel (abc) A {abc) 



Gf c) (Q 2 , Sl)S2 , S3 ) 



R-7T-K+ 

K s ir~K s 
K s -k~K l 
K-it°K° 



- cos 9 C 

cos 8 C 

- cos 6 C 

2 

COS^c 

2 

- cos 9 C 



Tf(g 2 )( v / 2 - 1) [v^r w (a 2 ) +rg( 5l ) 

Tf(g 2 )(v / 2 - 1) [v^T w (a 2 ) +4 1 ](5 1 ) 
Tf (g 2 )(V2-l)[T«( Sl )-T«( S3 y 
Tf(Q 2 )( v / 2 - 1) [2v^L,(aa) +2®(«s) 



TrVif" sin0 r 



if 7T 7T~ 



7T 



sin# r 



K°7i° V2sm6 c 



T { kI{Q 2 ) T$l( 8l )-T$l(8 2 ) 



(1), 



T$(g 2 ) [2rW( S2 ) + r$( ai ) + t$(* 3 ) 



,(1), 
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TABLE III. Predictions for the normalized decay widths T(abc)/T e and the branching ratios 
B(abc) for the various channels with (3k* = —0.135 in Eq. (JToj) . The contribution from the vector 
current is listed in column 3 and available experimental data are listed in column 5. The later are 
taken from [T^| , except for Kstt~Ks, where we quote the value for KshKs given in [19]. 



channel (abc) ^ 


T(abc)^ 


(pred.) 


V ^ e / anomaly 


{B{abc)f pred - ] [%} 


(B(abc)) iexpL) [%] 


K S 7i-K L 
K~ir°K° 


0.011 

0.011 

0.0027 

0.0058 

0.0090 


0.0045 
0.0045 
0.0008 
0.0029 
0.0032 


0.20 
0.20 
0.048 
0.10 
0.16 


0.20 ±0.07 
0.021 ±0.006 
0.12 ±0.04 


K~7T~7T + 


0.0080 
0.043 
0.054 


0.0007 
0.0043 
0.0058 


0.14 
0.77 
0.96 


0.09 ±0.03 
0.40 ±0.09 
0.41 ±0.07 


TABLE IV Same 


as Tab. 


ID for the Cabibbo suppressed decays with £ = 


— 33 in Ea (&\) 


channel (abc) ^ 


T(abc) j 


(pred.) 


(T{abc)\ {pred - ] 

V ^ e / anomaly 


(Biabc))^ [%} 


(B(abc)) {expL) [%} 


K~7T~7T + 


0.017 
0.077 
0.087 


0.0007 
0.0043 
0.0058 


0.30 
1.37 
1.55 


0.09 ±0.03 
0.40 ±0.09 
0.41 ±0.07 
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TABLE V. Same as Tab. 



Ill 



but with Pk* = 



-0.11 in Eq. @. 



channel (abc) 




/r(a&c)\ (precL) 

V ^ e J anomaly 


{B{abc)) {pre±) [%) 


(B(abc)f xpL) [%] 


K--n~K + 

K S TT'K S 

K S 7T-K L 
K-n°K° 


0.0106 
0.0106 
0.0026 
0.0055 
0.0085 


0.0042 
0.0042 
0.0007 
0.0027 
0.0030 


0.19 
0.19 
0.046 
0.098 
0.15 


0.20 ±0.07 
0.021 ±0.006 
0.12 ±0.04 


7r n°K- 


0.0076 
0.041 
0.052 


0.0007 
0.0041 
0.0056 


0.14 
0.74 
0.93 


0.09 ±0.03 
0.40 ±0.09 
0.41 ±0.07 


TABLE VI. Same 


as Tab. [II but with j3 K * = in Eq. (10). 


channel (abc) 


(T{abc)\ {pred - ] 

v r e J 


(T{abc)\ (pred - ] 

V ^ e / anomaly 


(£(a6c)) (precL) [%] 


(B(abc)) iexpL) [%} 


K~TT~K + 
Ks7T~Ks 

K s n-K L 
K-n°K° 


0.0084 
0.0084 
0.0020 
0.0044 
0.0067 


0.0034 
0.0034 
0.0006 
0.0022 
0.0024 


0.15 
0.15 
0.036 
0.078 
0.12 


0.20 ±0.07 
0.021 ±0.006 
0.12 ±0.04 




0.0060 
0.035 
0.045 


0.0005 
0.0035 
0.0048 


0.11 
0.62 
0.81 


0.09 ±0.03 
0.40 ±0.09 
0.41 ±0.07 
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